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We explore the salient features of the 'Kitaev ladder', a two-legged ladder version of the spin-1/2 
Kitaev model on a honeycomb lattice, by mapping it to a one-dimensional fermionic p-wave super- 
conducting system. We examine the connections between spin phases and topologically non-trivial 
phases of non-interacting fermionic systems, demonstrating the equivalence between the spontaneous 
breaking of global Z2 symmetry in spin systems and the existence of isolated Majorana modes. In the 
Kitaev ladder, we investigate topological properties of the system in different sectors characterized 
f"*) ' by the presence or absence of a vortex in each plaquette of the ladder. We show that vortex patterns 

\ can yield a rich parameter space for tuning into topologically non-trivial phases. We introduce a 

new topological invariant which explicitly determines the presence of zero energy Majorana modes 
at the boundaries of such phases. Finally, we discuss dynamic quenching between topologically 
non-trivial phases in the Kitaev ladder, and in particular, the post-quench dynamics governed by 
tuning through a quantum critical point. 
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I. INTRODUCTION 



It has been recognized recently that non-interacting fermionic lattice systems can exhibit a rich variety of phases 
^ I based on their topological properties. On the other hand, it has been long known that a class of spin systems can be 
mapped to non-interacting fermionic systems. Here, we exploit the mapping to tailor the recent insights in fermionic 
systems to the context of spin systems, in particular, to a system which we dub the 'Kitaev ladder'. The Kitaev ladder 
is a quasi one-dimensional (ID) analog of the Kitaev model on the honeycomb lattice a model which has been 
extensively studied for its topological aspects [l-Q; the ladder system shares many of the important features of the 
^ ] parent two-dimensional (2D) system. The ladder has two legs with spin-1/2 degrees of freedom at each lattice site, 
O • highly asymmetric couplings between nearest-neighbor spins, and conserved Z2 degrees of freedom on each plaquette 
O ^ [ associated with the presence or absence of a vortex. We will show that this ladder can be mapped to the celebrated 
one-dimensional p-wave superconducting system. 

In pioneering work [loj . Kitaev showed that a ID p-wave superconducting system can exhibit topologically dis- 
^ , tinct phases characterized by a topological Zo index. This index characterizes the topology of any one-dimensional 
■ particle-hole symmetric topological insulator [ll| . Topologically non-trivial phases are characterized by the presence 
of localized zero energy Majorana modes at the ends of an infinitely long system having open boundary conditions; 
5l5 ' kinetic energy and superconducting pairing, so to speak, conspire to split the Dirac fermions comprising the system 
. \ into their 'real' and 'imaginary' parts. Here, we discuss the manner in which these specific topological aspects and 
. phases translate to the context of spin chains. Topologically non-trivial phases in the fermionic system may be associ- 
ated with symmetry-broken phases in the spin systems, and the trivial phases with paramagnetic phases. As observed 
in [1, 3, we find that the presence of isolated Majorana modes at the ends of a long system with open boundary 
conditions is intimately related to spontaneous breaking of a global Z2 symmetry. We build on this observation by 
explicitly demonstrating the connection between spontaneous broken symmetry in a spin system and Majorana modes 
in the corresponding fermion system. While the behavior of spin expectation values and correlations is generally re- 
lated to those of fermions in a complex, non-local fashion, the language of topology provides an immediate and novel 
Ci • perspective on spin phases. 

The connection between topologically non-trivial and magnetic phases is generic to a large class of spin systems. 
What distinguishes the Kitaev ladder is the set of Z2 degrees of freedom associated with different vortex sectors. 
In the parent Kitaev honeycomb system, while most studies have focused on the ground state sector corresponding 
to the absence of a vortex in every plaquette, studies of other vortex sectors have been sparse. Here, the ladder 
system provides a simple prototype for studying the physics of difi^erent vortex sectors. The distribution of vortices 
on plaquettes can, in the corresponding fermionic system, be encoded in the sign of the chemical potential on each 
lattice site. We explore the phases realized in a range of periodic vortex sectors and find novel conditions for the 
existence of topologically non-trivial states. For instance, unlike in the regular ID p-wave superconducting system, 
some sectors require a critical amount of superconductivity (reflected by the magnitude of the superconducting gap 
parameter) to host non-trivial states. 

Of late, motivated by applications to topological quantum computational schemes [l^ , there has been a surg e of 
interest in seeking out systems and geometries that can realize and manipulate isolated Majorana modes [3, [TMl [13- 



l3l| . Some analogous studies in the Kitaev honeycomb system have identified Majorana modes bound to vortices and 
schemes for their manipulation [isj . In the ladder system, the isolated modes, as opposed to being present at vortices, 
completely parallel the ID p-wave superconducting system in being present at the interface between topologically 
trivial and non-trivial segments. Moreover, we find that the periodic patterns mentioned above provide a new route 
for finding phases and configurations that support these modes. In principle, for translationally invariant systems, such 
phases can be characterized by a Z2 topological index that considers the Berry phase accumulated by the eigenvectors 
in traversing the first Brillouin zone. In practice, we find that for complex periodic patterns, such a treatment proves 
to be rather involved, and can be replaced by the evaluation of a much more direct topological invariant derived from 
the equations of motion. Our method generalizes the identification of zero energy plane waves by Wen and Zee [11] to 
the case of evanescent modes. We use this scheme to pinpoint several different configurations of periodic patterns that 
yield localized Majorana modes in the bulk or at the ends of the Kitaev ladder. Our analysis also provides an alternate 
route for isolating Majorana modes in ID p-wave superconductors by applying appropriate periodic potentials. 

The manipulation of the Majorana modes requires dynamically changing the parameters associated with the Hamil- 
tonian describing the system. We consider dynamics from two angles. Based on the parallel between spin and fermionic 
systems, we briefly outline the first steps necessary for realizing schemes in the Kitaev ladder that are analogs of those 
recently proposed in the ID p-wave superconductor. Our second study of dynamics entails quenching between topo- 
logically trivial and non-trivial regions, an operation that is necessary for several dynamic schemes. The quench 
dynamics is interesting in its own right from the perspective of non-equilibrium quantum critical phenomena, and we 
find that it results in the generation of residual energy or defects whose density scales as 1/ y'r in the thermodynamic 
limit, where l/r is the quench rate. This defect production, we believe, also has implications for finite sized regions 
and could contribute to quantum decoherence. 

Our presentation is arranged as follows. In section |lll we recapitulate the key features of the ID p-wave supercon- 
ductor and its topological aspects. We then show its connections to spin physics in the context of the extensively 
spin-1/2 XY chain subject to a transverse magnetic field. In section Hill we introduce the model of interest, the 
Kitaev ladder, map it to the ID p-wave superconductor and outline techniques for isolating Majorana modes in this 
system. In section IIV| we present various vortex sectors in the Kitaev ladder system, the phases realized in these 
sectors, and the conditions for isolating Majorana modes. We discuss dynamical aspects in section IVl and, in section 
IVII the manner in which our studies and findings would translate back to the ID p-wave superconductor. 



II. CONNECTIONS BETWEEN ID p-WAVE SUPERCONDUCTORS AND SPIN CHAINS 

A. The Id p-wave superconductor: recapitulation 

The ID p-wave superconducting system of spinless fermions explored by Kitaev [l^ is described by the tight-binding 
Hamiltonian 

^ = E + -fn+lfn) + A [Un+l + fl+jX ) " M {flfn - 1/2)] , (1) 

n 

where w is the nearest- neighbor hopping amplitude, A the superconducting gap function (assumed real), and /i 
the on-site chemical potential. The translationally invariant system can be diagonalized in the momentum basis, 
/j. — -j= fne-~^^^ and shown to have the particle-hole symmetric dispersion 



ujk = ± y(2wcosfc + ^j)2+4A2sin^fc. (2) 

As shown in Fig. [U the system hosts distinct phases depending on the parameter |/i|/2u;; the system is gapped save 
for the regions demarcating the phase boundaries. The weak {intra-siie) pairing phase |/i|/2w < 1 (phases I and II in 
Fig. [ij is topologically non-trivial while the strong pairing phase (phases III and IV) is topologically trivial. 



Topologically trivial versus non-trivial phases are distinguished by the absence versus presence of boundary Ma- 
jorana modes, which can be visualized in a simple fashion in the extreme strong and weak pairing limits as follows. 
Consider decomposing the Dirac /-fermion above in terms of two Majorana fermions /„ = (a„ + ifo„) /2, where the 
Majorana fermions respect the relations = 6^^ = 1, = {&„,6m} = 2(5„i^„, — 0. In terms of the 

Majorana fermions, the Hamiltonian in equation ([l} becomes 

= ^ ^ [{-w + A) a„6„+i + (u- A) 6„a„+i - ^la^hn] , (3) 

n 
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FIG. 1: Zero temperature phase diagram of the ID p-wave superconducting fermionic system as well as the spin-1/2 XY chain 
in a transverse magnetic field. Phases I and II denote topologically non-trivial/ferromagnetic phases while phases III and IV 
denote topological trivial/paramagnetic phases. Inset: Dispersion relations near various gap closures at phase boundaries. 




as represented in Fig. Ufa). 

In the extreme weak pairing limit A = w > 0, /i = 0, the tight-binding Hamiltonian in ([1} couples the Majorana 
fermions in the staggered fashion shown in Fig. [2jc). A change in sign of A corresponds to shifting the staggered 
pattern by a lattice site, thus exchanging the roles of the a and b fermions. The superconductivity induced anomalous 
term and the normal hopping term conspire to separate a Dirac fermion into its Majorana components, leaving an 
isolated Majorana mode at each end. In this limit, equation ^ decouples into the form '^n=ii'^ + ^)ifnfn ~ 
where /„ = (a,i +i5„-|-i)/2 corresponds to Dirac fermions composed of the pairs of linked Majorana fermions depicted 
in Fig. HJc). This then yields two Majorana fermions h\ and cat which arc isolated at the two ends and which do not 
appear in the Hamiltonian since these modes have zero energy (loj . In the extreme strong pairing limit A = w = 0, 
/i 7^ 0, the Majorana fermions a„ and &„ arc pairwise connected as shown in Fig. [^b); thus no Majorana fermions 
are isolated. 

The presence of an energy gap in each phase ensures that slightly changing the couplings from these extreme limits 
does not alter these topological aspects. In the weak pairing phase, it is thus still possible to define appropriate linear 
combinations Ql = X^'^n'^n Qb,~ Pnbn (with a„ and (3,-1 real) that are Majorana modes bound to the two 
ends of the system and that become isolated in the thermodynamic limit. 



B. Mapping spin chains to a p-wave superconductor 

The concepts discussed above can be investigated in the context of spin chai n phys ics. As an example, consider the 
extensively studied spin-1/2 XY chain subject to a transverse magnetic field [33l435| . The Hamiltonian is given by 

N-l N 

H = - + Jv«+i) - (4) 

n—l n—1 

where a" denote the Pauli matrices. The spin-1/2 operators can be mapped to Majorana fermions using the Jordan- 
Wigner transformation [H, a„ = Jlj^i '^j '^n^^n = 11^=1 '^j for 2 < n < A^, ai = erf and hi ^ a\. The 
resultant Hamiltonian exactly maps to the ID p-wave superconductor Majorana Hamiltonian of equation Q with 
the identification w ^ -\- Jy^ A -f-)- Jj, — J^,, /z o 2h, and with the interchange of a„ and 6„ on even sites. Hence, 
Fig. [T] also represents the phase diagram for the spin system. Phases I and II are ferromagnetic and have non-zero 
ground state expectation values of and , respectively. Phases HI and IV arc paramagnetic with both and (t|J 
having zero expectation values in the ground state. The extreme limits shown in Fig. [2] can be understood in terms 
of spin physics. The case of Fig. ^c) corresponds to the values of the couplings Jx 0, Jy = 0,h = and thus to an 
Ising-type ferromagnetic (or antiferromagnetic) ground state along the x-axis. The case of Fig. [2Kb) corresponds to 
Jx = 0, Jy = 0,h and thus to spins polarized along the z direction. 



3 



(a) a, *2 ^3 a, 




FIG. 2; Pictorial representation of the couplings between Majorana fermions described by the Hamiltonian in equation Q 
for the cases of (a) general couplings, (b) the strong (intra-site) pairing limit A = w = 0, fj,^0 which does not possess 
Majorana modes (topologically trivial), and (c) the weak pairing limit A — —w, /i = which hosts Majorana modes at the 
ends (topologically non-trivial). 



C. Spontaneous symmetry breaking and isolated Majorana modes 

As discussed above, the topologically non-trivial phases in the p-wave superconductor, associated with isolated 
boundary Majorana fermions, map onto the ferromagnetic phases in the spin chain. We argue here as follows that the 
presence of these isolated fermions corresponds to a Z2 symmetry in the spin phases which gets broken spontaneously. 
As is obvious in the extreme limit of Fig. [5Jc) , the energetics in phases I and II provides — 1 constraints for a system 
of A^ spins. Hence, the ground state is constrained to two degenerate states reflecting a Z2 symmetry. Returning 
to the language of fermions, the degeneracy is associated with the isolated modes Ql = J2 cx-nO-n and Qb = ^ Pnbn 
which do not participate in the energetics. Explicitly, the Dirac fermion formed by the two boundary Majorana 
fermions, ^ = (Ql + iQR)/2 can have occupation number or 1. Thus, picking one of the corresponding degenerate 
states |0)^( or (or any linear combination u|0)* -I- u|l)>if) amounts to breaking the Z2 symmetry of the ground 
state and picking one of the two energetically available spin configurations. In phases III and IV, however, energetics 
poses A^ constraints and the system is confined to a unique ground state. 

To explore the spin physics in the symmetry broken phase in terms of the isolated Majorana modes, consider the spin- 
1/2 algebra formed by {Ql, Qr, —iQLQn}, where iQlQr = 24'^^!' — 1. The two eigenstates of any of these operators 
(those of the —iQlQr being |0)* and are orthogonal to one another and are each equally valid configurations 

in the ground state. The particular symmetry broken choice immediately determines the expectation values of these 
operators and their corresponding Majorana and spin operators. Specifically, if we consider the two eigenstates of Qi, 
IV'±)lj we can use the relationship {(3,a„} = 2q!„ to show that {'4!+\an\'4'+) l = otn and (V'- |a„| V-)-L = —OLn for all n. 
Furthermore, since {Q, 6™} = 0, we see that {'>p±\hn\'4'±) l = for all n. The physics behind these statements can be 
made more transparent, for instance, in the limiting case Jx > \ Jy\ and h = Q. Here, the system is in the Ising-type 
phase whose possible degenerate ground states correspond to spins that are either primarily aligned or anti-aligned 
along the x direction. In fact, \iP±)l are exactly these two states. Given that is identically ai, we can easily 
evaluate its expectation value in one of the states to get (^/i+jcrf = -y/l — [Jy/J^Y, as ought to be true for the 

system pointing primarily along x. 

The non-local nature of the Majorana modes is consistent with the spontaneously broken symmetry being global. 
Energetics locks all spins into one of two possible configurations where the Jordan- Wigner string connects the chain 
in a non-local, highly entangled fashion. A particular choice of states in the basis of the isolated Majorana modes 
thus corresponds to picking one of these highly entangled states. It should be noted that a knowledge of which one 
of these states is chosen, for instance gives limited information in that expectation values can be evaluated 

only for non-local spin operators corresponding to Majorana modes a„ and &„, except for one or two operators, such 
as af above. This knowledge does provide evidence for long-range order. For example, our result for (V'-i-luf 
in the example above can be compared to that of Lieb, Schultz and Mattis [s^. It is shown in Ref. [s^ that, for an 
open chain with A^ sites (where A^ is even and 00) in a case which is effectively that of the Ising limit above, the 
end-to-end two-point correlation functions are given by (-0+ 1 erf crj^ 17/1+) = 1 — {Jy / J^Y and {il]j^\a\a^j^\il]^) = 0. In the 
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limit iV — >• 00, we expect clustering to hold, so that | (^/i+l erf |?/'+)| = \{'ilj+\afa^\ip+)\^^'^ ■ We therefore find complete 
agreement with the results in [s^. We note that our derivation is much simpler than the one in Ref. [s^, purely 
making use of the end Majorana modes, and our results can be easily generalized to take into account the transverse 
magnetic field h in Eq. ^ which was not considered in [s^. Thus, analysis of the isolated Majorana modes can be 
a valuable means for gaining insight into spin chain physics. 



III. THE KITAEV LADDER SYSTEM 



A. The model 



Equipped with the connection between a ID p-wave superconductor and spin chains, we turn to the object of 
our attention, the Kitaev ladder system. This system consists of a single (modified) strip of Kitaev's original two- 
dimensional honeycomb system. The model is thus ladder-like and each plaquctte is a square instead of a hexagon. 
As shown in Fig. |31 the two-legged ladder hosts spin-1/2 degrees of freedom at sites (n, u) and (n, I), where u, I denote 
the upper and lower legs, respectively. The Hamiltonian is given by 

(W-l)/2 

^ — ^ [Jx {<^2n-l,u'^2n,u + '^2n,l'^2n+l,l) 
n=l 

N 

+ Jy (o-L,«<^L-n,« + 4n-i,/4n,;)] + J2 J^<,u<p (5) 

11=1 

where cr^ denote the Pauli matrices respecting the usual commutation rules, 



2n-l,u 2n,u 2n+\,u 




2n-\J 2n,l 2n+\,l 2n-\,l 2n,l 



FIG. 3: The Kitaev ladder, the associated couplings Jx, Jy and Jz, and the vortex operators, W, 



As in the 2D honeycomb lattice system, each plaquette has an associated operator Wn which commutes with the 

Hamiltonian and is of the form <yn,u(^n+i,u(^^n+i,i^^n+i,i " ^^^^ (^n,u'^ n+i, u'^n i'^n+i,i ^ Since = 1, 
the eigenvalues of Wn are ±1; hence these invariants provide a set of Z2 quantum numbers characterizing different 
sectors of the Hamiltonian. These operators, commonly referred to as vortex operators, reflect the absence/presence 
{Wn = ±1) of a vortex in each plaquette. The richness of this system, its phases and its topological properties derive 
from the extensive space of the Z2 configurations. 



B. Jordan- Wigner transformation 

The Kitaev ladder/honeycomb system is endowed with a special set of couplings which enable the Jordan- Wigner 
transformation to provide a local fermionic Hamiltonian 0, 0] • As in the previous section, this transformation gives a 
direct connection between the Kitaev ladder and the ID p-wave superconductor. The mapping relates spin operators 
to Majorana fermionic operators: 

a-n = Snj/uO-n.i/u^ ^" ^ ^^"^ " odd/cvcu, (6) 

Cn = Sn,l/u(^nJ/u^ ^" = ^^n,u/l(^n^u/V ^^^^ " odd/cveu. 

Here, assuming that each leg contains TV sites with open (rather than periodic) boundary conditions, the string 
operator Sn.u = Ilm^i '^m u runs from the left to the right along the top leg in Fig. |3]to the (n — l)-th site, while 
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Snd = Ylrn=i u Y\m=n+i '^m I ^^^^ complctcly along thc top leg left to right and then back along the bottom leg 
to the (n + l)-th site. 

In terms of these Majorana operators, the Hamiltonian is of the form H = 
'^^^{—iJxa„bn+i + iJybnCin+i — JzO-nbnC-ndn)- The role of the vortices is encoded in the Jz term. This can 
be seen by noting that Wn = — (ic„(i„)(ic„+i(i„+i) and defining s„ = ic„d„ = Yim^n^"^'^^'^^ ■ Then, the 
Hamiltonian, as depicted pictorially in Fig. [31 takes the form 

-ff = j ^ i-JxC-nhn+l + JyKan+l + JzSnanbn) , (J) 
n 

where s„ = ±1 depending on the configuration of vortices, namely Wn = — s„s„+i. Now, to make a connection with 
the p-wave system, one can once more define Dirac fcrmions /„ = (a„ + ib„)/2 as in the XY spin chain. The vorticity 
can be captured by another set of Dirac fermions, gn = (c„ + idn) /2; s„, which we coin 'site-polarity', reflects the 
occupancy of these fermions via s„ = 2gl^g^ — 1 and determines the sign of the chemical potential for the /-fermions 
In terms of the Dirac fermions, the Kitaev ladder also maps onto the ID p-wave superconductor form, 

n 

- A (/„/„+! + fi^ji ) - /i„ (/t/„ - 1/2)], (8) 

where the identification 

w^Jx + Jy, A^Jy-Jx, and ^„ f> 2s„ (9) 

has been made. Thus, thc /-fermions, as in the XY spin chain, participate in thc dynamics, while the (7-fermions 
encode the vortex configurations. Each set of fermions spans a Hilbert space of size 2^, together comprising the 
Hilbert space of size 2^^ thus accounting for all the degrees of freedom of the original 2N spins of the ladder system. 

Compared to the XY spin chain of the previous section, one main difference in the mapping to the p-wave super- 
conductor is that the sign of the chemical potential is position dependent. Here, we restrict our studies to periodic 
patterns. For s„ with period p, we can employ the decomposition 

p-i 

s„ = ^ (^s^e'^'^'^/P + s;e-'2'^«"/Pj . (10) 

q=Q 

In the momentum basis {fk = -^^n /n^ ''^"), the Dirac Hamiltonian describing the Kitaev ladder then takes the 
form (up to an overall constant) 

0<fe<71- 

p-1 

0<fc<27r q=0 

The other difference between thc ladder and the XY spin chain is in thc degrees of freedom. Here, the ladder 
Hamiltonian provides A'^ — 1 constraints on the energetics for a ladder of 2N spins. An additional — 1 constraints 
arc given by the configuration of vortices (note that another Z2 symmetry is present because the vortex eigenvalues 
of Wn are invariant under a global sign change of the s„'s). These constraints together connect all 2 A'' spins. Hence, 
once again the topologically non-trivial phases, namely those bearing isolated boundary Majorana modes, can be 
associated with spontaneous symmetry breaking in thc spin language; a detailed study of the spin physics, as initiated 
for the XY spin chain, is in order. 

C. Techniques for identifying isolated Majorana modes 

The fermionic Hamiltonian in Eq. ([8]) provides a starting point for exploring various conditions under which bulk 
or boundary Majorana modes may appear. In this section, we develop a formalism which enables us to identify these 
modes in a range of vortex sectors. We avail ourselves of the transfer matrix method which has been used extensively 
in ID systems [37l . [ssj and is well suited to study bound states. The presence of bound Majorana states is governed 
by the growth or decay of the eigenfunctions of the transfer matrix. As a way of discerning the existence of these 
states, we employ an invariant which is constructed by analytically continuing plane wave states. 
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1. Transfer matrix approach 



The transfer matrix can be constructed by employing the Heisenberg representation and deriving the equations of 
motion from the Hamiltonian of Eq. ([7|) for the time dependent Majorana modes a„ = a„e~*'^* and 6„ = /3„e^'"*. 
These equations are of the form 

[w - A) a„_i + («; + A) q;„+i - /x„a„ = -iujfin, 
- (w + A)/3„_i - (w - A) /3„+i + M«/3« = -jwttn, (12) 

where we have invoked the identification in Eq. ([9]). 

To identify Majorana modes, we focus on cj = corresponding to real values of a„ and The equations for a„ 
and Pn decouple, and they can be written in the transfer matrix form 

'^"+1 \ - A ( \ „.lnorp A - ( sq ' 



O-r, 



= A„ ( 1 , where A„ = ( A+^« A+» j _ (13) 



Similar expressions hold for the /3„'s since the respective transfer matrices arc related by _B„ — A^^. Knowing the 
behavior of the a- modes thus completely determines that of the 6- modes. 

This transfer matrix formulation enables us to study the growth versus decay of modes at the boundary of a 
finite piece of the ladder or at the interface between two parts of the ladder in different phases. In particular, for a 
homogeneous system having all couplings and s„ constant, this behavior is determined by the eigenvalues of any An, 
while for a region having periodicity p in the s„'s, it is determined by the eigenvalues of the matrix 

Ap = Ap A.p-i ■■■ A2 Ai. (14) 

Majorana modes bound to the ends of a finite chain require that both eigenvalues of be either smaller or greater 
than unity in magnitude. The former case corresponds to an a-modc localized at the left end and a 6-modc at the 
right end. The conditions on the eigenvalues ensure that the localized modes can simultaneously respect constraints 
at the boundary and normalizability. For instance, if the chain goes from n = 1 to 00, the condition to be imposed 
on the Majorana mode at the left end is that ao = and ai = c, where c is a non-zero constant to be fixed by 

the normalization. The corresponding vector ^ ^ = ) ^'^ written as a linear superposition of the two 

eigenvectors of ^p. We then see from Eq. (fT3|) that q;„ ^ as ?i — > 00 if both the eigenvalues of .4p are smaller than 
unity. 

Bound Majorana states in the bulk of the chain can be realized by juxtaposing two regions corresponding to 
different phases. The highlighting feature of this system is its ability to dial in to different phases by changing the 
vortex patterns. We thus consider the interface of two semi-infinite patterns P (n < 0) and Q {n > 0) having 
periodicities p and q, respectively, but with the same values of the couplings w, A, and fi. Considerations similar 
to those above enable us to identify two situations in which there exists exactly one normalizable Majorana mode 
for n tending to both 00 and —00. These correspond to (i) Ap having both eigenvalues larger than 1 in magnitude 
or both smaller than 1 and Aq having exactly one eigenvalue less than 1 in magnitude, or (ii) Ap having exactly 
one eigenvalue larger than 1 in magnitude and Aq having both eigenvalues less than 1 or both smaller than 1 in 
magnitude. Other situations yield either no bound mode or a pair of concurrent Majorana modes. 



2. Topological invariant 

The relevant features of the transfer matrix structure can be gleaned by a simple consideration of its analytic 
properties. This analysis extends the application of the topological invariant introduced by Wen and Zee [3^ to the 
study of evanescent modes. Our construction of the invariant is derived from f{z), the characteristic polynomial of 
Ap. Given that Ap is a real 2x2 matrix, a single root inside the unit disk must lie on the real axis. Thus the index 

= - sgn(/(l)/(-l)), where f (z) ^ dct{Ap - I z) , (15) 

is equal to 1 if and only if ^p has exactly one eigenvalue with a magnitude less than 1. This case of = 1 indicates 
a topologically trivial state, i.e., one with no end Majorana modes. However, v = —1 indicates the existence of end 
Majorana modes reflected by Ap having both eigenvalues less than or both greater than 1 in magnitude. The marginal 
case ;/ = would imply a closure of the bulk gap. In fact, v can only change sign if a bulk gap closes. 
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The invariant v is topological in that it is derived by invoking the continuity of f[z) to infer its behavior inside 
the unit disk based on its values on the unit circle. We believe that this index, derived using physically transparent 
methods, can be reduced to Kitaev's invariant based on the Pfaffian when generalized to our ladder model. The 
invariant v can be related to that of Wen and Zee [s^l by noting the form v = (— 1)"^"'"^, where Uj is the number of 
zeros of /(z) for which \z\ < 1; in Ref. [l^, the equivalence Uf = ^ §^^^^1 dz f'{z)/f{z) was used to identify plane 
wave zero modes. Furthermore, f{z) was replaced by the dispersion which can also be done in our case but makes 
the analysis more involved. Finally, the invariant v can be related to the topological invariant of [ll| in that they 
both locate zero modes but the former bypasses the need to identify the eigenvectors of Eq. (|lip and the Berry phase 
across the Brillouin zone. 

IV. VORTEX SECTORS AND ISOLATED MAJORANA MODES 

We employ the transfer matrix approach and the invariant u to show that a diverse set of patterns in the site- 
polarities Sn, or equivalently, vortex sectors, can host isolated Majorana modes. We focus on three classes of vortex 
sectors (i) the full vortex sector defined by s„ = 1, (ii) the zero vortex sector defined by s„ = (— 1)", and which is 
associated with the ground state sector in the parent 2D model, and (iii) patterns of s„ having higher periodicity. 

A. Full vortex sector: s„ = 1 

The full vortex sector maps directly to the XY spin chain in the transverse field as well as the homogeneous ID 
p-wave superconductor. Isolated Majorana modes have been well-studied in the latter case and are expected to exist at 
the ends of wires in the topologically non-trivial phases I and II in Fig. [T] and at the boundary between a topologically 
non-trivial phase and topologically trivial phase in the phase diagram of Fig. [1] We expand on these arguments to 
illustrate the techniques developed here and to provide a comprehensive description of the Kitaev ladder. In the next 
subsection, we also show that this sector corresponds to the ground state sector in a certain parameter range. 

1. End modes 

The condition for Majorana modes to exist at the ends of the Kitaev ladder is given by ly = —1, where ly is the 
topological invariant of equation (jl5p . For the full vortex sector, we have a uniform site-polarity, s„ = 1 for all n; the 
corresponding site-independent transfer matrix is given by 



A=i^ A+.„ j ^ (16) 

which yields 

f{z) = det(^ -Iz) = z'-( -J^) z ( ^] . (17) 



wj \A 
Hence, 

= - sgn(/(l)/(-l)) = sgn(|^|-2|^|), (18) 

for A -I- w 7^ 0. A consideration of the location of the zeros of f{z) shows that an a (b) Majorana mode is localized 
to the right (left) side of the system if 2w > and A > 0. If 2w < — 1^| the positions of the a and b Majoranas are 
switched. Switching the sign of A also interchanges the Majorana modes (see Fig. 

2. Interface modes 

As an illustrative example, consider the bound modes obtained by juxtaposing two regions in phases I and III on the 
left and right, respectively, in the phase diagram of Fig. [TJ The bound state structure at the interface can be derived 
by considering the closing of the gap at the phase boundary. Close to this region, the dispersion in equation ^ takes 
the form Wfc = ±2y/ A^{dk)^ + at fc « tt, where m{x) — /i/2 — w. Here, the continuum version of the equations 
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of motion, equation (|T2|) . is given by dxa{x) = m{x)a{x) and dxb{x) = —m{x)h{x). In this case, the Majorana mode 
that satisfies normalizability takes the form b{x) cx e~ '^^ ). In fact, this solution corresponds to that of the 
celebrated Jackiw and Rebbi model [s^ as applied to Majorana fermions. The instances of a single bound Majorana 
mode at an interface are exhausted by the cases involving the juxtaposition of phases I/II < 2\w\) and phases 
III/IV (|/i| > 2|w|). In terms of the invariant form in Eq. p7)) . these cases correspond to juxtaposing a topologically 
non-trivial (i' = — 1) phase with a topologically trivial phase (i^ = 1). 




+ + + + + + +- + - + - 




v=-l 



FIG. 4: Examples of isolated Majorana modes for the Kitaev ladder system. Generically, Majorana modes are hosted at (a) 
the ends of ladders with nontrivial topology [v = —1) or at (b) the interface between regions with v — —1 and v = 1. In the 
case of the uniform sector Sn — 1 with 2|ii;| > |/^|, individual Majorana modes are hosted at (c) each end of a ladder. Other 

examples of end Majorana modes include (d) the alternating sector s„ = ( — 1)" with 2|A| > and (e) the sector P — +H 

with 2^y\wA\ > \fi\. 



B. Vortex-free sector: s„ = ( — 1)" 

In the 2D honeycomb system, the vortcx-free sector is special in that it corresponds to the ground state sector for 
all values of spin couplings. This condition, however, holds only for a range of parameters in the Kitaev ladder. To 
elucidate this, consider the limit of Jz large compared to Jx and Jy. In the language of a p-wave superconductor, 
let us assume that ^ is large and positive. We then find that the perturbative change in the ground state energy 
compared to the w ~ A = limit is given by the second order expression 

AEo = - Y K(l-SnS»+i) + A2(1 + s„s„+i)] . (19) 

n 

If > A^, the ground state energy is lowest if s„s„+i = —1 for all n, i.e., if s„ = (—1)"; this is the sector with no 
vortices. If < A^, however, the ground state energy is lowest if s„s„-|_i = 1 for all n, i.e., if s„ = ±1 for all n; 
this is the full vortex sector. Another difference is that unlike in the 2D system, reduced dimensionality renders the 
Kitaev ladder system to be gapped everywhere except along phase boundaries. 

For the vortex free sector, we can identify the range of parameter space which yields topologically non-trivial phases. 
Towards this end, we identify the period 2 transfer matrix 

fj. A — w \ / A— m \ 

A+w A+w 1 ( A+iii A+tu 1 (or)\ 

1 J \ 1 J ' ^ ' 
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and the characteristic polynomial 



/(z)=det(^-/z) 



z + 



A- 



(w + A)" ) \^ 
Our topological invariant (for w + A 7^ 0) takes the form 

V ^ - sgn(/(l)/(-l)) = sgn(|/i|-2|A|) 



(21) 



(22) 



In this case, the topological nature of the phase depends on the relative magnitude of the superconducting order 
parameter and the chemical potential, as opposed to the hopping integral and the chemical potential as in the 
uniform case, cf. Eq. ([TS]) and Eq. (|22p . That the point separating the topologically trivial and non-trivial regions 
derived from = is equivalent to gap closure can be gleaned from the dispersion for s„ = (—1)", 



= 2 



2 



[li? - A^)cos(2fc) ± A^tsinfc 



(23) 



Specifically, consistent with Eq. ((22)) . the gap closes at fc = ±7r/2 for 2|A| = Once again, as in the full vortex 
sector, the isolated Majorana modes at the end of a finite chain or at the interface of two regions can be identified 
using Eq. [22] (see Fig. |4j;) . Thus we see that the imposition of a periodic pattern can give rise to qualitatively different 
physics in that here the topological nature of the system depends on the magnitude of A. 



C. Higher period sectors 



The consideration of periodic patterns in the site-polarity s„'s provides a zoo of configurations for realizing topo- 
logically non-trivial phases and interface Majorana modes. We can once again invoke the topological invariant v for 
any periodic structure P characterized by the transfer matrix „4_p of Eq. ([T4l) . In general, we have 



(24) 



i/p = - Sgn(detylp + 1 + mTrylp) 



m— ± 



sgn(|Trytp| - |dety^p + l|). 



It follows from these expressions that any two chains which are cyclic permutations of each other must have the same 
topology, as we would expect physically. It is also true that the transformation s„ — > — Sn for all n will not change the 
topology of the chain. Although giving a full analytic expression is not practical, it may be noted that the detailed 
features of a pattern of the site polarities s„ that control its topology enter Eq. ([24]) through the quantity Tr^p. Some 
insight can therefore be gained by an expansion in certain limits. For example, if < 2\w ± A|, a useful expansion 
of TrAp in powers of ^j? /{w"^ — A^) is (for p even) 



w 



p/2 



2(u.2 - A2 



p/2 

E 

L m— 1 



'p/2 



(25) 



Thus, in this case the quantity ^X]m=i •^2™^ (^J2n=i ■^2n-i^ represents the most relevant feature of P in determining 
its topology. As a specific example, taking Ap ~ A1A2A3A4 for a generic period 4 pattern (51,32,53,54), we get an 
expression for which is a function of (51 -I- 53)(52 -I- 54), 51525354, and the couplings constants w, A, and fi. Thus 
from the standpoint of topologically non-trivial phases, there are two distinct period 4 patterns, namely -I- -1- H — and 
+ H . For instance, the application of equation to the pattern P = + -\ yields 



vp = sgn (/.j'^ — 4|wA|) 



(26) 



as illustrated in Fig. |4f. Similar constraints are also obtained for the other patterns, demonstrating that patterns of 
higher periodicity can have phase boundaries involving more intricate dependencies on the values of fi, w, and A. 
Thus, we have shown that the topology and existence of isolated Majorana modes can be controlled by configurations 
of vortices in the Kitaev spin ladder. 
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V. DYNAMICS 



Having charted out the topologically non-trivial phases and associated Majorana physics in the Kitaev ladder, we 
move on to dynamically tuning between topologically trivial and non-trivial regions. We briefly visit this issue in the 
context of manipulating isolated Majorana modes and then present a comprehensive study on dynamic quenching 
between these regions. 



From the point of view of topological quantum computation, a system hosting zero energy Majorana modes can 
in principle store non-local quantum information. In 2D systems, gate operations on the degenerate subspace can be 
achieved by the dynamic interchange of such modes, exploiting the non-Abelian nature of such exchanges. Recent 
proposals have explored schemes for transporting isolated Majorana fcrmions in ID p-wave superconductor wires and 
also braiding them in a network of quantum wires including T-junctions (l6l. [l9j. These protocols are based on the 
tunability of the bulk topology and the controlled growth or contraction of regions of non-trivial topology. 

The Kitaev ladder is well-suited to such a scheme given the many different ways of controlling the topology of the 
system. Primarily, moving an isolated Majorana mode would amount to dynamically moving the interface between 
topologically trivial and non-trivial regions, such as those shown in Fig. |4l One method of doing so would be to 
dynamically change the spatial profile of the couplings Jx^ Jy, and Jz which determine the phases in the different 
regions, for instance, as shown in the phase diagram of Fig. [1] Alternatively, Majorana modes may be manipulated by 
changing the vortex patterns which, as discussed in section IIVI determine the topology of a region. The presence or 
absence of a vortex is determined by the site-polarity s„, or, equivalently, the f/-fermion occupation number discussed 
in section llll Bl These occupation numbers can be controlled by a combination of operators which do not appear in the 
Kitaev ladder Hamiltonian in equation (0). Examples include Ji{<J2n,u'^2n+i,u-<^2n,i<^2n+i,i) = 2«^i(ffn5„+i 

and J2{cr^n,u^2n+i,u + o-fn,iO'2«+ij) = '^^Mglgn+i + ffnsi+i)- A possiblc mcans of reversing the site-polarities on 
adjacent sites would be by the application a "tt pulse" involving either of these operators. Schemes for exchanging 
Majorana modes can be developed using such transport and, in principle, given that the system is a ladder, as opposed 
to a wire, the need for T-junctions can be eliminated. 

Translating between the language of qubits hosted by Majorana modes and spin physics is a rather subtle task 
requiring further study. As discussed in section III CI a direct correspondence exists between isolated Majorana 
operators and certain spin operators. Care needs to be taken to ensure that the relevant degrees of freedom arc not 
affected by local fields, and hence, environmental effects. Naively, this seems possible given that the Jordan- Wigner 
string is non-local, and at the same time, the final read out could be more accessible than in fermionic systems by the 
use of spin-sensitive probes such as magnetic cantilevers. Thus, spin systems may afford an interesting laboratory to 
explore topological phases of matter and the ideas of topological quantum computation. 



The issue of quenching a topological ID system through different quantum phases is particularly germane for schemes 
which involve moving bound Majorana modes localized by shifting topologically non-trivial segments [l6l [l9j. The 
dynamics associated with a slow quench across a quantum critical point (QCP) has generated considerable recent 
interest in its own right j40l - l52| . although its application to topological systems is as yet sparse. To recapitulate some 
of the basic features, quenching typically involves initializing a system in the ground state of a given Hamiltonian 
and then varying a parameter of the Hamiltonian at a finite rate 1/r so as to take it across a QCP. The prominent 
aspect of such a quench is that one ends up with a finite density of defects, or a residual energy (defined below), 
which scales as a power of l/r. In d dimensions, the defect density or residual energy obeys the Kibble-Zurck scaling 
form 1/t'^'^/'^^'^~^^\ where and z are respectively the correlation length and dynamical critical exponents of the QCP 
[40l - |45j . These exponents are defined by the divergence of the correlation length as |A — Ad"'' and the relaxation time 
as |A — Xc\~^'^ , respectively, as a parameter A in the Hamiltonian approaches a critical value Ac. 

Here, we consider quenching for the case extensively discussed in previous sections, namely, the full vortex sector 
(s„ = 1) in the Kitaev ladder which maps directly onto the ID p-wave superconductor. In this case, in terms of the 
Majorana fermions, in the momentum basis, theHamiltonian ([5]) is of the form 



A. Majorana manipulation 



B. Quench dynamics 




(27) 
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where fl denotes the momentum vector (/^ f-k)- As detailed earlier, the system respects the dispersion ujk = ± 



(2wcosfc + /i)^ +4A2sin^A: and the resultant phases are as shown in Fig. [T] We study two kinds of quenches: 
(A) between the two possible topologically non-trivial phases / and // in Fig. [l] and (B) between a topologically 
non-trivial and topologically trivial phase, such as between phases / and ///. We assume a linear quench in which one 
of the parameters of the Hamiltonian, to be specified below, varies in time as a(t) = wot/r, so as to take the system 
through the relevant QCP. We quantify the effect of the quench using the residual energy per site, Er, attained at 
the final time, defined as 

t,N^oo N\a{t)\ ^ ' 

Here, N denotes the number of sites, (if) / denotes the expectation value of the Hamiltonian in the final state reached, 
and i?o is the ground state energy at t — >■ oo. We will study the dependence of Er on t for r ^ 1; we expect this to 
vanish in the adiabatic limit r — > oo. 

For both kinds of quenches, as shown in the dispersions in Fig. [1] the system consists of two bands which are 
symmetric about the Fermi energy. In the ground state, the band below the Fermi energy is occupied and the one 
above is unoccupied. Well within one of the phases, the system is gapped. The quench takes the system through a 
gapless point at a QCP (as described below) and then into another gapped phase. For any given set of fc-modes in 
Eq. (P7)) . the quench couples the two states in the subsystem respecting /^t/^, + f^i^f_i^ = 1. Thus, the relevant sector 
contributing to the quench dynamics comprises of the states having occupation numbers |rifc, n_fc) equal to |1, 0) and 
|0, 1). In this basis, the quenches for cases (A) and (B) take the form 

TiAkit) = {2wocosk + ^lo) t'' + a(t)(2sinfc) T^, 

Hskit) = {2wocosk + a{t)) T~ + Ao(2sinfc) r"^, (29) 

where wq-, fJ-o and Aq correspond to some fixed parameter values, and r° denote the Pauli matrices. Thus, case (A) 
involves tuning through A = for fixed |/io/(2wo)| < 1 and case (B) through fj, = 2wo for fixed Aq/wq. These take 
the system through a gapless point which lies at wave vectors kc — zLcos~^[—fio/{2wQ)] in case (A) and k^ = and tt 
in case (B). In both cases, the critical exponents are given by ^ = = 1, since the gap is (i) linear in |fc — /ed at the 
QCP a = 0, and (ii) linear in \a\ when a is close to zero, for k = kc- 

For the linear quench a(t) = wot/r, th e dy namics in both cases maps exactly to the well-known Landau-Zener 
problem, as with a host of other quenches [4l| - |48l . [53j. In particular, for case (B), the probability pk of ending in an 
excited state at time t and the net residual energy, for which each sub-system contributes Apk , are given by 

dk 

Pk = exp[-7r(2Aosinfc)V/u;o], and Er = / — Apk- (30) 

Jo 27r 

The probability pk is largest for the low-energy modes near fee = and tt where the gap vanishes at the QCP. In 
the limit r ^ wq/Aq, the residual energy is dominated by these modes; expanding around k^ gives a Gaussian 
integral and the power-law form Er ^ l/r^/^ which is consistent with the Kibble-Zurek power dv /{zv -\- 1) given that 
d ^ z = V = 1. The same power- law is obtained for case (A), although the expression for the excitation probability 
Pk looks a bit different from case (B). 

In principle, quenching can be studied in the other vortex sectors discussed above. For most sectors, the problem 
becomes analytically intractable even for a case in which A = and only /i is varied in time linearly. As shown 
by two of us in previous work (s^ . an interesting situation arises when a pattern of s„ is endowed with a higher 
symmetry in that some matrix elements of the Hamiltonian in equation (|lip are equal to zero. For instance, in the 

sector in which the signs of s„ form the period 4 pattern (+ H ), we find that the low-energy modes at fc = and 

TT (which dominate the generation of residual energy) are coupled via a single intermediate high-energy state lying 
at fc = — 7r/2 and 7r/2. Using second order perturbation theory to eliminate the high-energy states, we find that the 
effective time-dependent Hamiltonian governing the low-energy states is given by 

^-ff - \ 2k -Ct^/T^ ) ' (^^) 

where the variable fc has been redefined to lie close to 0, and c is a constant. Using a scaling analysis, one can then 
show that the residual energy goes as Er ^ [U- Further work would involve studying such anomalous scaling 

for quenches into topologically non-trivial phases involving A 7^ 0. 

More generally, if the most relevant states participating in the quench are not directly coupled, but instead are 
coupled via q intermediate steps, it can be argued that the residual energy, for a linear quench, scales as Er 
\ I / ii+'i) ^ Thus, the quench dynamics in these sectors is also controlled by the QCP, resulting in a power-law 
scaling of the relevant quantities with respect to the quench rate. 
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VI. APPLICATIONS TO THE ID p-WAVE SUPERCONDUCTOR 



Given that most of our analysis of the Kitaev ladder involves mapping to the ID p-wave superconductor, several of 
our results can be directly applied to the latter system and related ones. 

(i) Periodic potentials and topologically non-trivial phases:- One of our main results is that vortex patterns in the 
Kitaev model can alter the topology of the system. The cases that we considered involved periodic patterns in the 
signs s„. The primary effect of a periodic pattern is to introduce more bands in momentum space, an effect which 
can equally well be achieved by the presence of a superlattice or a periodic externally applied potential. Thus, we 
predict that in the p-wave superconducting system, an alternate means of changing the topological nature of a wire 
and the associated Majorana end modes is the application a periodic potential. We have developed a straightforward 
technique for identifying these features even for complex dispersions based on a topological invariant directly derived 
from equations of motion. Our analysis provides detailed results for the topologically non-trivial phases and isolated 
Majorana modes in the special case of the external potential controlling only the sign of the local chemical potential 
on each site. Our results show that such a potential can alter the physics in qualitative ways, for instance, in the case 
of an alternating site potential, enabling the superconducting gap function to act as a tuning parameter for achieving 
topologically non-trivial phases. 

(ii) Quench dynamics:- The quench analysis above provides a word of caution for schemes which entail dynamically 
changing the topological nature of wire segments. In the thermodynamic limit, our analysis shows that no matter 
how slowly the system is tuned from one phase to another, the presence of the QCP results in defects whose density 
has a power-law dependence on the tuning rate l/r. Thus, except in the extreme adiabatic limit, the quench 
is bound to generate defects which could cause quantum dccohcrcnce. Even in the case of changing a finite segment 
from a topologically non-trivial to a topologically trivial region, it should be kept in mind that the eigenstates of the 
initial system span its entire length and a local change, particularly involving a gapless point, could affect the entire 
system and produce a finite density of defects. In fact, the route to equilibrium upon quenching a finite region is in 
itself an interesting problem which could place stringent constraints on proposed dynamical schemes. 

VII. CONCLUSIONS 

In conclusion, we have introduced and analyzed a ladder version of Kitaev's honeycomb model. By employing a 
mapping to a ID p-wave superconducting system, we have borrowed from the insights provided by this system in 
terms of topological features. We find an intimate connection between spontaneous symmetry breaking and isolated 
Majorana modes; this connection provides us with an alternative view of topological order and offers an attractive 
avenue for the investigation of a wide class of spin systems. 

We have performed a detailed study of the ladder system in the presence of vortex arrays and shown that the presence 
of vortices can dramatically alter the topological properties of the system. We have identified various conditions and 
configurations of vortices that result in the isolation of Majorana modes, states which are currently being considered 
as prime candidates for the building blocks in topological quantum computation. We have instigated the study of 
quench dynamics across quantum critical points in these topological systems, discussing the power-law scaling of the 
residual energy density; we believe that this avenue unveils important physics from the perspective of non-equilibrium 
quantum critical phenomena as well as in realistic treatments of the dynamic manipulation of topological entities in 
quantum computational schemes. 

Our findings are germane to various physical systems that arc currently under study for their topological properties. 
Several of our results are directly applic able to quantum wires which, by various methods, can effectively host proximity 
induced superconducting order [l6l4l9l [23l [25l - [29l [Slj. In the past years, several proposals have focused on realizing 
the Kitaev honeycomb system in cold atoms trapped in an optical lattice jssl [Sq ; our results can be immediately 
translated to these systems. Furthermore, our studies show that it is possible to build a dictionary between topological 
aspects of non-interacting fermionic systems and ordering in spin systems; an extensive dictionary would enable us to 
apply insights found in one physical system to the other. 
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